NASA Technical Memorandum 104367 


//&>?/ 

(?i l \ 

Analysis of Some Compliance 
Calibration Data for Chevron-Notch 
Bar and Rod Specimens 

(NAS A- TM- 104367) ANALYSIS OF SOME N91-22594 

COMPLIANCE CALIBRATION DATA FOR 
CHEVRON-NOTCH BAR AND ROD SPECIMENS (NASA) 

14 p CSCL 20K Unci as 

G3/39 0011691 


Thomas W. Orange, Raymond T. Bubsey, William S. Pierce, 
and John L. Shannon, Jr. 

Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

Symposium on Chevron-Notch Fracture Test Experience 
sponsored by the American Society for Testing and Materials 
Indianapolis, Indiana, May 6, 1991 





ANALYSIS OF SOME COMPLIANCE CALIBRATION DATA 
FOR CHEVRON-NOTCH BAR AND ROD SPECIMENS 


Thomas W. Orange, Raymond T. Bubsey,* William S. Pierce,* and John L. Shannon, Jr. 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

This paper presents a set of equations describing certain fracture mechanics 
parameters for chevron-notch bar and rod specimens. They are developed by fitting 
earlier compliance calibration data. The difficulty in determining the minimum stress 
intensity coefficient and the critical crack length is discussed. 


INTRODUCTION 

This paper presents a set of equations describing certain fracture mechanics 
parameters for chevron-notch bar and rod specimens. They are developed by fitting 
previously reported experimental compliance calibration data. Their use will facilitate 
the testing and analysis of both brittle metals and the tougher ceramics. The 
equations present the various parameters in forms suitable for determining fracture 
toughness from maximum load, for determining the crack-extension resistance curve 
(R-curve), and for setting instrument sensitivities. The data encompass the entire 
range of the specimen geometries most commonly used. 

We first discuss briefly the background of the chevron-notch specimens and the 
experimental data to be used. Then we present a more extensive discussion on some 
particular characteristics of the chevron-notch specimens and their practical 
application. The fitted equations are presented and their fitting accuracies are 
discussed. Finally, problems in determining the minimum stress intensity coefficient 
and the critical crack length are discussed. 


SYMBOLS 

a Crack length (measured from load line) 

a 0 Distance from load line to tip of chevron 

a m Crack length at which Y* is minimum 

B Specimen thickness 

C Specimen compliance, C = EBV/P 

C' Compliance derivative, dC/da 

D Diameter (rod), D = B 


'Retired. 


E Elastic (Young’s) modulus 

K, Opening-mode stress intensity factor 

Plane-strain fracture toughness for chevron-notch specimens 
P Applied load 

V Crack mouth opening displacement 

W Width 

Y Dimensionless stress intensity factor for a crack In a chevron notch, 
Kpw' /Z /P 

Y* m Minimum value of Y as a function of a 

a a/W 

«o 

a m/ W 


BACKGROUND 

The chevron-notch specimens are fairly recent additions to the field of fracture 
mechanics. Consequently they do not have the same historical background of 
extensive stress intensity and displacement analysis as do the more common 
specimen types. But, like the earliest specimen types, we can develop useful 
expressions using experimental compliance data. 

Compliance data for the chevron-notch bar [1] and rod [2] specimens were 
previously reported. In each paper, one fitted equation was presented relating the 
minimum stress intensity factor to the initial crack length and to the specimen 
dimensions. A later paper [3] reported additional data for specimens having smaller 
initial crack lengths and also revised the previous equations to cover the wider range 
of crack lengths. But those equations alone are not sufficient for all analyses and 
tests involving high-toughness ceramics. To make them more complete and useful, a 
new set of generalized equations are presented in this paper. These equations are 
developed by fitting curves to the existing data. They are usable over a wide range of 
specimen dimensions. 


CHARACTERISTICS OF CHEVRON-NOTCH SPECIMENS 

For most common fracture test specimens, the dimensionless stress intensity 
factor (Y) increases continually with increasing relative crack length (a/W). But due to 
the wedge shape of the un-notched material in the chevron-notch specimen, the 
corresponding factor (V*) reaches a minimum, denoted Y m , as the crack length 
reaches a value denoted a m . The values of Y m and a m are functions of specimen 
dimensions and notch geometry only and are independent of material properties. 

If the material being tested has a crack growth resistance curve which increases 
rapidly to a relatively constant plateau (known as a "flat" R-curve), instability will occur 
at a=a m and P=P max . Then the fracture toughness (KJ can be calculated from 
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and no other test measurements are necessary. 

For some materials (even some ceramics), however, the R-curve does not reach a 
plateau but continues to increase with increasing crack extension (a "rising" R-curve). 
For such materials eqn. (1) does not apply and it may be desirable to determine the 
complete R-curve. In this case ASTM Test Method E 561 [4] may be used for 
guidance. If crack mouth opening displacement (CMOD) is measured during the test 
(as in E561) and appropriate compliance relations are available, one can calculate the 
instantaneous crack length. From crack length and toad, one can calculate the crack 
extension resistance as 


k x = Y* 




( 2 ) 


A plot of crack extension resistance against crack advance is the R-curve. 


PROCEDURE 

Experimental 

The experimental procedure is described in detail in Refs [1,2]. The complete data 
are presented in Ref [5]. At least three replicate tests for each crack length were 
averaged to obtain the data reported here. For each specimen, 7 to 15 crack lengths 
(depending on the initial crack length) were tested. 


Basic data reduction 


Analysis of the data is based on the following equation [1] 
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and its derivative with respect to a. In [1-3] the logarithms of the basic compliance 
data (C = EBV/P) were fit with a fourth degree polynomial in a. The fitted curve was 
differentiated and the values of V calculated from Eq 3. 

In Ref [1] the reported values of / m and a m corresponded to the minimum of that 
fitted curve. In Ref [2], / m and a m were determined in the same way, but the data 
range was restricted to seven points symmetrical about the value of a m found by the 
first fitting. Ref [3] used still another procedure. Seven points were selected by the 
previous criterion. Then a fourth degree polynomial was fit to the logarithms of the 
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compliance derivatives. That second polynomial was used to calculate / m and a m . 

In the process of verifying these calculations, some general concerns arose 
concerning procedures for determining y* m and a m . These will be discussed later. 


Development of generalized equations 

The following expressions are useful for computing the plane strain fracture 
toughness when the material has a relatively "flat" R-curve. 

a* = A,, + + A 3 cto < 4 ) 

and 

Ym * + + Bjflto (5) 

These were developed by first fitting third-degree polynomials in a 0 for each specimen 
type (bar or rod) and each value of W/B. Then the coefficients of the intermediate 
polynomials were in turn fit to a second-degree polynomial in W/B to produce the final 
forms of Eqns (4) and (5). Values of the coefficients for Eqns (4) and (5) are given in 
Tables 1 and 2, respectively. 

An expression for determining the relative crack length a as a function of 
measured displacements is 

(6) 

a = C 0 + C X U + CzU 2 + C 3 U* + C A U* 


where U is the Saxena and Hudak form [6] 


u = 


l 


| EBV j 1 / 2 


+ 1 


The coefficients for Eqn (6) are given in Table 3. This equation lends itself to 
computer-controlled fracture toughness testing since the subcritical crack growth can 
be determined from automated load and displacement data acquisition. 

When the relative crack length a is known, the stress intensity factor / and the 
dimensionless compliance EBV/P can be computed from the following expressions: 

Y* u e A> + Ol« ♦ ♦ Dj * 1 ♦ D t m* ( 7 ) 


and 
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Bo ♦ By* + By* 2 ♦ By* 3 * K 4 «* 


The coefficients for Eqns (7) and (8) are given in Tables 4 and 5, respectively. 


DISCUSSION 


Generalized equations 

Eqn (4) fits the calculated values of o m within 0.01 3W for the bar specimens and 
within 0.006W for the rod specimens. Eqn (5) fits the calculated values of within 
1.0% for the bar specimens and within 2.7% for the rod specimens. 

Within the ranges of a and a 0 specified in Tables 3 to 5, Eqn (6) fits the measured 
values of a within 0.003W for the bar specimen and within 0.002W for the rod 
specimen; Eqn (7) fits the calculated values of / within 2.9% for the bar specimen and 
within 2.1% for the rod specimen; and Eqn (8) fits the measured values of EBV/P 
within 1 .4% for both the bar and the rod specimen. 

Table 3 of Ref [7] gives values of V* m and a critical slope ratio r c . That ratio is the 
ratio of the compliances corresponding to a m and a 0 . For specimens with W/B =2.0, 
the values of Y* m computed from Eq (5) for both the bar and rod specimens are within 
0.6% of those in Ref [7]. The critical slope ratio computed from Eq (8) is within 1% for 
the bar specimen but is 7.8% low for the rod specimen. 


Problems in determining and a m 

The method of data analysis used in Ref [3] was not given explicitly and could not 
be determined directly from archival records. In attempting to verify the numerical 
analysis (by duplication), several methods were tried. Each produced a significantly 
different value for a m for the same data set, and this is a problem that should be 
discussed. 

The problem is inherent in the chevron-notch specimen. It is due to the same 
characteristic that makes it desirable, namely the fact that / has a minimum. For 
example, assume that we have a function f such that 

EBV/P = f(a) 

where f includes the data transform (if any) and a fitting function. Substituting this into 
the derivative of Eq (3) and eliminating non-zero terms we have 

o = i — f'{ a n ) - (9) 

a m - a 0 


where V and f are the first and second derivatives and a m is the root of this equation. 
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Unlike simpler specimens, we need to determine the second derivative as well. This 
presents a strong challenge to the analyst. 

Fig 2, from Ref [1], shows the typical variation of V* with a for different values of 
o 0 . Experimental compliance data would be expected to scatter about these lines. It 
is apparent from this figure that for a short initial crack (say, a 0 =0.2) Y* m will be 
relatively insensitive to the method of curve fitting but a m will be very sensitive. 
However, for a long initial crack (say, a 0 =0.5) the opposite will be true. 

Thus if the primary objective of the test is to determine K*. the initial crack length 
should be short. This is the case in Ref [7]. However, a long initial crack length is 
preferable if the critical crack length is important for, say, fractographic purposes. 

It should be pointed out that numerical a nalyse s (i.e., finite element or boundary 
integral methods) are subject to the same problem, although to a lesser degree. 
Discrete pairs of (/, ag) for several initial crack lengths must be fitted with a function 
to calculate a minimum. Three pairs are required, more would be preferred. 


CONCLUSIONS 

The equations presented here are in forms suitable for several purposes in fracture 
testing with chevron-notch specimens. They encompass the range of specimen 
geometries most commonly used and provide a good fit with the basic compliance 
data. The inherent difficulty in determining the critical crack length from compliance 
measurements is discussed. 
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Table l--Coefficients for Eq 1 


Specimen Coeff. 

Expression 

Ao 

Bar A, 

a 2 

A 3 

- 0.110 +0.354(W/B) -0.088(W/B) 2 
0.268 +1.628(W/B) -0.400(W/B) 2 
1.637 -6.358(W/B) +1.872(W/B) 2 
0.075 +4.462(W/B) -1.508(W/B} 2 

Ao 

Rod A, 

a 2 

A 3 

0.147 +0.089(W/B) -0.026(W/B) 2 
0.358 +1 . 1 50 (W/B) -0.096(W/B) 2 
2.860 -5. 190(W/B) +0.770(W/B) 2 
-3.610 +5. 100(W/B) -0.800 (W/B) 2 

Range: 1 .5<(W/B)<2.0, 

0<ap<0.5 



Table 

2- -Coefficients for Eq 5 

Specimen 

Coeff. 

Expression 


Bo 

-17.03 +29. 94 (W/B) -5.0(W/B) 

Bar 

B, 

-116.00 +141.60(W/B) -29.6(W/B) 


b 2 

1131.00 -1304. 00 (W/B) +342.0(W/B) 



-1351.00 +1654.00(W/B) -443.2(W/B) 


Bo 

5.47 + 6. 29 (W/B) +2.46(W/B) 2 

Rod 

Bi 

-65.93 + 72. 62 (W/B) -5.62(W/B) 2 


b 2 

622.00 -659. 80 (W/B) +146.10(W/B) 2 


B 3 

-541.40 +629.10(W/B) -135.20(W/B) 2 


Range: 1 .5<(W/B)<2.0, O^a^O.S 


Table 3--Coeff1cients for Eq 6 


Specimen 

Coeff . 

Expression 





Co 

3.09 

-24.12a,, 

+57 

•12a 0 2 

Bar 

Ci 

-109.30 

+1 227.00O0 

-2 876 

.OOa 0 2 

W/B=l .5 

C 2 

1 908.00 

-22 216.00oq 

+51 286 

. OOa 0 z 


C 3 

-14 900.00 +168 580.00O0 - 

381 240 

. 00a o 2 


C4 

41 390.00 - 

451 059.00oq +987 080 

1 . 00a o 2 


Co 

2.08 

-8.74op 

+16. 

93ao 2 

Bar 

Ci 

-63.31 

+540.00a 0 

-1 019. 

OOa 0 2 

W/B=2.0 

C 2 

1 086.00 - 

11 296.00oq +20 043. 

00a o 2 


C3 

-9 327.00 +98 493.000^ -158 690. 

00a* 


C4 

28 430.00 -284 970.00a 0 +366 330. 

OOa 0 2 


Co 

0.672 

+4.85a 0 


-23.93a 0 2 

Rod 

Ci 

25.670 

-361.90oo 

+1 

624.00a o 2 

W/B=1.5 

C2 

-858.000 

+9 512.00a o 

-39 

580 . 00a 0 2 


C3 

9 219.000 

-105 260.00O0 

+411 

440 . 00a 0 2 


C4 

-35 145.000 

+417 050.000o 

-1 550 

300.00a o 2 


Co 

0.896 

+7 . 24a 0 


-26.5a 0 2 

Rod 

Ci 

21.800 

-590.40a 0 

+2 

087. Oa 0 2 

W/B=2.0 

C2 

-1 192.000 

+17 166.00a o 

-58 

980. Oa 0 2 


C3 

16 772.000 

-213 330.00O0 

+713 

640. Oa 0 2 


c« 

-78 837.000 

+961 870.00oq 

-3 146 

400. Oa 0 2 


Range: 0. 18<aQ<0.22, ao<a<0.8 



Table 4--Coefficients for Eq 7 


Specimen 

Coeff. 

Expression 


Do 

3.329 +1 .026oq +78. 210^ 

Bar 

D, 

-0.812 -58.080O0 -334. 400^ 

W/B=l .5 

d 2 

-2.061 +265.260oo +461. 40a,, 2 


d 3 

4.350 -417.12000 -156. lOa^ 


d 4 

0.349 +219.800O0 -65.55oo 2 


Do 

4.308 -4.757a 0 +83. 770^ 

Bar 

D, 

-6.529 -19. 190a 0 -358. 70a,, 2 

W/B=2.0 

d 2 

16.630 +172.000a o +483. lOc^ 2 


d 3 

-22.170 -313.000a,, -151. lOa^ 


d 4 

13.220 +173.70000 -72. 710^ 


Do 

-2.28 +106. 3a 0 -567. Oa 0 2 +1 062a 0 3 

Rod 

D, 

29.61 -582. 5a 0 +3 100. Oa^ -5 830a 0 3 

W/B-1.5 

d 2 

-60.17 +1 167. Oa 0 -6 070. Oa 0 2 +11 589a 0 3 


d 3 

52.60 -1 022. Oag +5 051. Oa 0 2 -9 869a 0 3 


d 4 

-15.00 +337. 8a 0 -1 517. 0a/ +3 042a 0 3 


Do 

2.19 +41.95a 0 -263. 7oo 2 +749. Oa 0 3 

Rod 

D, 

-1.41 -128.20a 0 +945. 6a 0 2 -3 532. Oa^ 

W/B=2.0 

d 2 

22.30 +38.71a 0 -691. 5a* +5 784. Oa 0 3 


d 3 

-40.95 +179. 10a o -703.00^ -3 646. Oa 0 3 


d 4 

22.92 -131. 30oq +741. 2a 0 2 +607. 4a 0 3 


Range: 0.1<ao<0.35 (bar), 0 . l<a o ^0 . 40 (rod), a o <oc<0.8 


Table 5--Coefficients for Eq 8 


Specimen 

Coeff . 

Expression 


E 0 

2.850 -6.4800 +61.56a„ 2 

Bar 

Ei 

1.177 +26.59oo -349. 30a,, 2 

W/B=l .5 

e 2 

9.650 -8.37a,, +708.90oo 2 


e 3 

-16.240 -62.60a o -597.00oo 2 


E 4 

10.450 +56 . 82a 0 +167.90oo 2 


E 0 

3.885 -17.750o +94.97oo 2 

Bar 

Ei 

-5.160 +123.20oo -624.20a 0 z 

W/B=2.0 

e 2 

34.270 -324.50oo +1 562.00a o 2 


e 3 

-52.330 +386.80oo -1 756. 000^ 


E 4 

27.950 -173.20a 0 +741.40a 0 2 


E 0 

3.91 -23. 18a„ +138. 4a 0 2 -91.84a 0 3 

Rod 

Ei 

-10.01 +237.70a 0 -1 356. OaQ 2 +1 325.00a 0 3 

W/B=l .5 

e 2 

51.60 -758.80a 0 +4 284. Oa 0 2 -4 777.00a 0 3 


e 3 

-74.66 +969.60a 0 -5 480. Oa 0 2 +6 516.60a 0 3 


E 4 

37.83 -433.00oo +2 464. 5a 0 2 -3 043.00a 0 3 


Eo 

2.92 +0.28oo +26.67a 0 2 +33.27a 0 3 

Rod 

Ei 

1.68 +28.52oo -336.00a 0 2 +111.40a 0 3 

W/B=2.0 

e 2 

20.59 -135.84a 0 +1 157.00a 0 2 -918.40a o 3 


e 3 

-39.16 +218. 10a 0 -1 581.00a 0 2 +1 581.00a 0 3 


E 4 

22.64 -115.78a 0 +756.60a 0 2 -837.20a 0 3 

Range: 0. 

l<Oo<0.35 (bar) 

, 0.1<Oo<0.40 (rod), Oo<a<0.8 
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